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ABSTRACT 

Using multi-wavelength imaging from the Wide Field Camera 3 on the Hubble Space 
Telescope we study the stellar cluster populations of two adjacent fields in the nearby 
face-on spiral galaxy, M83. The observations cover the galactic centre and reach out to 
~ 6 kpc, thereby spanning a large range of environmental conditions, ideal for testing 
empirical laws of cluster disruption. The clusters are selected by visual inspection to be 
centrally concentrated, symmetric, and resolved on the images. We find that a large 
fraction of objects detected by automated algorithms (e.g. SExtractor or Daofind) 
are not clusters, but rather are associations. These are likely to disperse into the 
field on timescales of tens of Myr due to their lower stellar densities and not due to 
gas expulsion (i.e. they were never gravitationally bound). We split the sample into 
two discrete fields (inner and outer regions of the galaxy) and search for evidence of 
environmentally dependent cluster disruption. Colour-colour diagrams of the clusters, 
when compared to simple stellar population models, already indicate that a much 
larger fraction of the clusters in the outer field are older by tens of Myr than in 
the inner field. This impression is quantified by estimating each cluster's properties 
(age, mass, and extinction) and comparing the age/mass distributions between the 
two fields. Our results are inconsistent with "universal" age and mass distributions of 
clusters, and instead show that the ambient environment strongly affects the observed 
populations. 
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1 INTRODUCTION 

The standard paradigm that has been built up over recent 
years is that most, if not all, stars form in clusters, and that 
the majority are disrupted due to the removal of the gas 
left over from the non-100% star formation efficiency (e.g. 
Lada & Lada 2003). This process, called 'infant mortality', 
is thought to be largely independent of cluster mass (Hills 
1980), and the timescale over which it operates is several 
internal crossing times of the embedded cluster. However, 
the importance of this infant mortality process depends cru- 
cially on how a "star cluster" is defined (Bressert et al. 2010, 
Gieles & Portegies Zwart 2011). 

Beyond this early phase, theoretical and numerical 



studies (e.g. Spitzer 1958; Henon 1961; Baumgardt & 
Makino 2003) suggest that the survival of a cluster should 
depend on its mass, with more massive clusters surviving 
longer, and environment, with strong tidal environments or 
high giant molecular cloud - GMC - densities causing higher 
disruption rates. However, observationally, the case is less 
clear, with reports of a strong dependence of disruption on 
cluster mass and environment (e.g. Lamers et al. 2005a, b) 
and no dependence on either parameter (e.g. Whitmore et 
al. 2007). Quantifying the effects of disruption in cluster 
populations is essential in order to 1) constrain the frac- 
tion of stars that are originally formed in clusters, which in 
turn may affect star formation theories (e.g., Bastian 2008; 
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Silva- Villa & Larsen 2011), and 2) use clusters to trace the 
star formation histories of full galaxies or individual regions 
within them (e.g., Pellerin et al. 2010). 

Lamers et al. (2005a) and Gieles, Lamers, & Porte- 
gies Zwart (2007) have suggested that cluster populations 
in the Galaxy and SMC have been shaped by cluster mass 
and environmentally dependent disruption, respectively (see 
Larsen (2006) for a thorough review of empirical mass de- 
pendent disruption laws). In this scenario, high mass clusters 
are more likely to survive to older ages, and environments 
with low GMC densities and/or small tidal fields are more 
conducive to cluster survival. Hence regions or galaxies that 
have a low GMC density and small tidal forces should con- 
tain a larger fraction of older clusters relative to regions with 
strong tidal forces. 

Based on observations of the Antennae galaxies (Fall 
et al. 2005,2009; Whitmore et al. 2007,2010), the SMC and 
LMC (Chandar et al. 2006; 2010a), and the inner regions of 
M83 (Chandar et al. 2010b), these authors have suggested 
that cluster disruption is independent of mass and environ- 
ment. They suggest that a certain fraction of clusters, Fmid, 
is disrupted every decade in age (with Fmid= 0.8 — 0.9). This 
strong cluster disruption means that the formation history 
of the population is a secondary effect, therefore popula- 
tions should follow "universal" age and mass distributions 
(although see Chandar et al. (2011) for specific examples 
where this does not seem to hold within M51). In this sce- 
nario, cluster disruption is driven largely by internal mech- 
anisms. 

Finally, Elmegreen & Hunter (2010) have suggested a 
time-dependent disruption scenario, where clusters are dis- 
rupted by the hierarchical interstellar medium from which 
they formed. This is similar to the Gieles et al. (2006) models 
where GMCs dominate the disruption process, however, in 
the Elmegreen & Hunter (2010) scenario, the GMC density 
experienced by the cluster is time-dependent (i.e. allowing 
a cluster to drift away from a gas rich environment), mean- 
ing that most of the clusters are destroyed at young ages 
(Kruijssen et al. 2011). This scenario predicts that cluster 
disruption may be independent of mass but strongly depen- 
dent on environment. 

Here, we attempt to address the environmental depen- 
dence (or lack thereof) of cluster disruption, using WFC3 
imaging of two adjacent fields within M83, a nearby (~ 
4.5 Mpc, Thim et al. 2003) face-on barred spiral galaxy. 
The central pointing covers the nucleus and inner spiral arm, 
while the outer field has an average galactocentric distance 
of ~ 2 times that of the inner field (2.5 kpc vs. 4.75 kpc, or 
0.25 vs. 0.47 R25 - Paturel et al. 2003). Due to its smooth (al- 
though warped) outer Hi distribution and extended UV disk 
(Thilker et al. 2005), it is reasonable to assume that M83 
has not undergone a recent strong interaction or merger, al- 
lowing us to assume a roughly constant SFR over the past 
0.5 — 1 Gyr, barring the inner nuclear region, which has 
an ongoing starburst (e.g. Harris et al. 2001 and references 
therein). The assumption of a constant star- formation rate 
over the past ~ 1 Gyr is the same as was made by Chan- 
dar et al. (2010b). Additionally, we have used the resolved 
stellar photometry and techniques presented in Silva- Villa 
& Larsen (2011) to derive the SFH in two fields at differ- 
ent galactocentric radii (centred at ~ 2.9 and 4.8 kpc from 
the galactic centre, however they are different fields than 



those used in the present work). We do not find an appre- 
ciable difference in the SFH over the past 100 Myr (after 
which the data become incomplete). This will be discussed 
in more detail in a forthcoming paper. 

Combining the two fields allows us to study a wide range 
of environmental conditions, perfect to study its effect on 
cluster disruption. Chandar et al. (2010b) have already used 
the inner field images to study cluster disruption in M83, 
and find evidence for mass and environmental independent 
disruption. Here we expand on their work by considering 
a broader range of environments, namely by studying the 
inner and outer regions of the galaxy simultaneously. Addi- 
tionally, we apply the same cluster selection techniques in 
both fields, which allows us to carry out a comparative study, 
with any effects of potential selection biases minimised. 

In a forthcoming paper (Bastian et al. in prep.) we will 
provide a detailed analysis of the dataset used here, along 
with a comparison with the catalogue of M83 clusters used 
by Chandar et al. (2010b), and a more detailed discussion 
on cluster disruption mechanisms. This is necessary in order 
to fully address the role of cluster mass in the disruption 
process. 



2 OBSERVATIONS 

2.1 Data and Cluster Selection 

For the present work we use Early Release Science data of 
M83 (GO 11360, PI O'Connell) taken with the WFC3 on 
HST of two adjacent fields in M83. The data for the inner 
field have been presented in Chandar et al. (2010b), and 
for the current study, we only make use of the observations 
in the F336W (U), F438W (B), F555W (V), F657N (Ha), 
F814W (I) filters. While we use the shorthand notation for 
the filters, we note that no transformations were carried out. 
The outer field data consist of imaging in the same filters 
(with the exception that the "V-band" was taken with the 
F547M filter), with similar exposure times, and will be pre- 
sented in more detail in a future work (Bastian et al. in 
prep). The data was taken fully reduced from the WFC3 
Early Release Science websitfl 

As a first step in creating a cluster catalogue we ran 
SExtractor (Bertin & Arnouts 1996) over the B, V and I 
images with settings chosen to select a large number of can- 
didates. These three catalogues were then cross-correlated, 
and ISHAPE (Larsen 1999) was run on the resulting list to 
search for resolved objects (FWHM > 0.2 pixels). The V- 
band image of each cluster candidate was then examined, 
and only resolved, centrally concentrated and symmetric 
sources were retained as clusters. Approximately 40-50% of 
the initial candidates were removed during this step. The 
vast majority (> 98%) of the objects remaining in our sam- 
ple after this step satisfy the definition of a cluster suggested 
by Gieles & Portegies Zwart (2011), namely that the age of 
the cluster is longer than a crossing time (based on radius 
and mass estimates). While this method of cluster selection 
naturally introduces a certain level of subjectivity, we found 
that it was necessary in order to remove unbound (i.e. fila- 
mentary) groups/associations and chance alignments of two 

1 http:/ /archive. stsci.edu/prepds/wfc3ers/m83datalist. html 
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or more bright stars. See Silva- Villa & Larsen (2011) for a 
more thorough discussion of selecting clusters. Details of the 
sample selection will be presented in a future work, however 
we note that the same selection procedure was used in both 
fields, so any differences should be physical and not due to 
any possible selection biasefl 

We then removed all clusters within 450 pixels (~ 
400 pc) of the galactic center, as the star formation rate 
is likely to have been highly variable there (e.g. Harris et 
al. 2001), and the detection limit is significantly worse than 
in the surroundings. Aperture photometry was then carried 
out on all the images with an aperture, inner and outer back- 
ground annuli or 5, 8, & 10 pixels, respectively. Photometric 
zeropoints were taken from the STScI website, with an ad- 
ditional 4% efficiency in the F336W filter (priv. comm. Ja- 
son Kalirai) taken into account due to better than expected 
efficiency. Aperture corrections for each filter were derived 
using a list of ~ 15 resolved clusters common to all filters. 
Comparison between our photometry and that of Whitmore 
et al. (2011, hereafter Wll) shows good agreement (devia- 
tions in colour < 0.03 mag, accounting for the F336W shift, 
which was not applied in Wll). 



2.2 Estimating Cluster Parameters and 
Comparison with Previous Works 

We estimated the age, mass, and extinction of each of the 
clusters by comparing the observed cluster magnitudes to 
simple stellar population (SSP) models. Two methods were 
used, in the first we adopted the GALEV SSP models (Ko- 
tulla et al. 2009) of 2.5 times solar metallicitj(f], and a 
Kroupa IMF and used the 3DEF fitting code (see Bastian 
et al. 2005). In the second method, we adopted the fitting 
procedure of Adamo et al. (2010a,b) and the Yggdrasil SSP 
models (Zackrisson et al. 2011), adopting the same stellar 
population parameters as the GALEV models. Both sets of 
models include nebular emission (e.g. Ha, [Oni]) which is 
helpful in distinguishing old clusters from young, highly ex- 
tincted clusters (e.g. Chandar et al. 2010b; Konstantopou- 
los et al. 2010) and contributes to the broad-band colours. 
Both methods gave consistent ages/masses for the clusters. 
In the present work we adopt the ages/mass derived using 
the Adamo et al. (2010a) methorfl 

We limit our analysis, when cluster parameters are fit, 
to those clusters that have masses in excess of 5 x 1O 3 M 
in order to minimise stochastic sampling effects that can 
severely affect the derived age/mass distributions (e.g. Maiz 
Apellaniz 2009; Silva- Villa & Larsen 2011). After applying 
this limit (and avoiding the inner ~ 450 pc of the galaxy) 
our final sample contains 381 and 370 clusters in the inner 
and outer fields, respectively. 

2 The full catalogues of observed and estimated cluster properties 
can be found on the CDS website. 

3 We adopt 2.5 times solar metallicity for both fields, as the 
metallicity gradient in M83 is such that we only expect a change 
in metallicity of A(12 + log(0 / H)) 0.2 in the mean of the two 
fields (Kewley et al. 2010). 

4 In a future work (Adamo et al. in prep) we will compare the 
results in more detail, and investigate systematic effects based 
on the combination of filters and models used (e.g. Adamo et 
al. 2010b; Reines et al. 2010). 
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Figure 1. A comparison of the estimated ages of clusters in the 
Chandar et al. (2010b) sample with those presented in this paper 
(adopting the ages from the Adamo et al. (2010a) method). Note 
that only clusters appearing in the Wll catalogue are shown. The 
dashed line shows a one-to-one correlation. The clusters marked 
with squares are those for which Wll assign older ages (their 
category 5a or higher) based on Ha morphology. For clusters 45, 
182, 232, and 294 our derived ages match those estimated by Wll. 



In Fig. \T\ we compare the ages we derived (using the 
Adamo models) vs. the modelled ages presented in Wll 
(taken from the catalogue of Chandar et al. 2010b). Overall, 
the agreement is extremely good. The three outlying clus- 
ters, where we assign significantly older ages than derived 
by Chandar et al. (2010b) all lack significant Ha emission, 
and their morphology suggests older ages (Wll). 



3 TRACING CLUSTER DISRUPTION 

To address the role of environment in cluster disruption we 
have used two methods, which are presented in this sec- 
tion in order of their dependence on models. The first, and 
least-model dependent, uses the observed colour distribu- 
tions directly. The second uses the ages and masses esti- 
mated through comparison with SSP models (see § 12. 2p to 
study how the age and mass distributions differ as a function 
of environment. 



3.1 Colour-Colour Diagrams 

In Fig. [2] we show the (U-B) vs. (V-I) colours of all clus- 
ters in the inner and outer fields. Note that the V-band 
corresponds to the F555W and F547M bands in the upper 
and lower panels, respectively. However, the y-axis uses the 
same combination of filters, and hence this is where we will 
concentrate our analysis. A clear difference between the two 
fields is the larger fraction of clusters in the outer field that 
have redder (U-B) colours, while the (V-I) colour distribu- 
tions are comparable. The significance of this result is easily 
seen when SSP models are overplotted on the data (red solid 
lines), which span from young (< 10 Myr) in the lower left 
to old (> 1 Gyr) ages in the upper right. Hence it is already 
clear from this plot that the average age of the clusters is 
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F547M-F814W 

Figure 2. Colour-colour diagrams for the inner (upper panel) 
and outer (lower panel) fields. Each cluster is shown as a filled 
circle, while the contours denote the number density of points at 
that colour. The (red) solid lines are the SSP model tracks used 
in Adamo et al. (2010b). For reference the models are marked 
with asterisks at 1, 10, 100 & 1000 Myr, from lower- left to upper- 
right. Note the differences in the distributions, with a significantly 
higher fraction of clusters in the outer field having colours consis- 
tent with being older. The average error in colour for our cluster 
sample is shown in the upper left of each panel, and the extinction 
vector is also shown. 

older in the outer field, suggesting that the cluster popula- 
tion this field is less affected by disruption. This conclusion 
is independent of any assumptions on cluster metallicity. 

In principle, this result could also be due to a higher ex- 
tinction in the outer field, but this is not a favoured scenario 
as 1) the inner regions host more gas/dust per unit surface 
area (e.g., Crosthwaite et al. 2002) and 2) we do not see 
as many reddened clusters (V — I > 0.8) in the outer field. 
In Fig. [5] we are showing all detected clusters in each field, 
although we note that mass and luminosity limited samples 
show the same distributions. 

3.2 Age and Mass Distributions 

Using the estimated ages and masses from § 12.21 we can di- 
rectly compare the age/mass distributions between the two 
fields. We chose to follow the method of de Grijs & Good- 
win (2008), also used by Chandar, Fall, & Whitmore (2010a) 
and Silva- Villa & Larsen (2011), namely to split each sample 
into three age bins, and plot the mass functions normalised 
to the linear age range of each bin. The results are shown in 
Fig. [3] Note that the mass functions extend to higher masses 
in the inner field. This is likely due either to size-of-sample 
effects (e.g. Hunter et al. 2003, Gieles & Bastian 2008) or 



through a higher upper cut-off in the mass function, i.e. M* 
if the mass function is described as a Schechter function 
(Larsen 2009, Gieles 2009). 

If there was no cluster disruption we would expect each 
of the three distributions in each panel to lie on top of 
each other. However, since the distributions are separated 
it shows that either the cluster-formation history has been 
increasing over the timespan shown (e.g. as may be expected 
in starburst galaxies, see Bastian et al. 2009), or more likely 
in this case, disruption has been acting on the population. If 
the differences between the distributions were due entirely 
to an increasing SFH over the past Gyr, an increase of a 
factor of 10 and 4 would be required in the inner, and outer 
fields, respectively. 

Additionally, we overplot the expected "universal" clus- 
ter distributions of Whitmore et al. (2007), Fall et al. (2009), 
and Chandar et al. (2010a,b) as dashed lines, normalised to 
the 10 — 100 Myr age bin. For this, we adopted a mass- 
independent fraction of clusters that are destroyed each 
decade in time, choosing Fmid= 0.9, and a cluster mass 
function of N(dM) oc M Q dM, with a = -2. We have 
chosen to normalise the "universal" distributions to the 
10 < age/Myr 100, in which case these distributions 
over and under-predict the younger (^ 10 Myr) and older 
(> 100 Myr) distributions, respectively. 

Comparing the upper and lower panels in Fig. [3] we see 
that the distributions of the three age bins are vertically 
spaced further apart in the inner field. The fact that the ob- 
servations lie nearly on-top of each other in the outer field 
suggests that disruption has had a relatively small effect 
on the cluster population far from the galactic centre. Ad- 
ditionally, we have carried out the same analysis for both 
fields, only using clusters with ages < 100 Myr (i.e. where 
the derived SFH based on the resolved stellar population is 
known to have been constant, see § 1), and come to the same 
conclusion. The observed differences between the two fields 
clearly shows the environmental dependence of cluster dis- 
ruption, with disruption affecting the inner field more than 
the outer field. 



4 DISCUSSION AND CONCLUSIONS 

Using two methods, we have shown that cluster disruption is 
dependent on the galactocentric distance (i.e. the ambient 
environment including GMCs and the tidal field) in M83. 
In the inner regions of the galaxy, where tidal effects are 
stronger and the density of GMCs is higher, we find evi- 
dence for higher levels of cluster disruption relative to the 
outer field. This result is based on a purely empirical com- 
parison of the cluster mass functions versus age in the two 
fields, and does not depend on any specific theoretical as- 
sumptions about the physical mechanism behind cluster dis- 
ruption. Our results could also be explained by differing clus- 
ter formation histories between the two fields, however the 
derived SFH based on the resolved stellar population does 
not support this scenario, at least for clusters that formed 
during the past 100 Myr. 

The observed differences between the inner and outer 
fields is in agreement with the mass and environmentally 
dependent disruption model of Lamers et al. (2005b), as 
the gas surface density in the inner region is higher by 
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Figure 3. The mass distributions for three age bins (normalised 
to the linear age range considered) for the inner (top panel) and 
outer (bottom panel) fields. The points represent the median mass 
in each mass bin, while the horizontal lines show the mass range 
covered by the bin. The dashed lines show the predicted "uni- 
versal" cluster distributions of the mass independent disruption 
scenario (Fmid = 0.9) normalised to fit the 10-100 Myr distri- 
butions. We have used bins with equal numbers of clusters (10 
per bin). The distributions in the upper and lower panels are 
clearly different, suggesting that disruption has been much less 
pronounced in the outer field than in the inner field. 



a factor of 3 — 4 than the outer region (Crosthwaite et 
al. 2002). The observed differences also agree, at least in 
a qualitative sense, with the mass-independent cluster dis- 
ruption model by Elmegreen & Hunter (2010) as well as the 
mass-dependent model presented by Gieles et al. (2006), as 
the ambient GMC density is expected in both models to 
strongly affect the disruption process. However, the "univer- 
sal" distribution (e.g. Whitmore et al. 2007; Fall et al. 2009; 
Chandar et al. 2010a,b), where mass and environmentally 
independent cluster disruption is the dominant factor in the 
population, clearly does not fit the data. 

In a future work, we will compare our catalogue of clus- 
ters in the inner field more thoroughly to that presented in 
Chandar et al. (2010b) in order to quantify the effects of 
cluster selection. This is necessary to fully address the role 
of mass in the cluster disruption process. Additionally, we 
will analyse the mass distributions of clusters and the results 
of the full three dimensional model fitting (4, Fmid, M*) to 
the observed distributions. 
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